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Abstract 

Simple derivation is presented of the four families of infinitely many shape invari- 
ant Hamiltonians corresponding to the exceptional Laguerre and Jacobi polynomials. 
Darboux-Crum transformations are applied to connect the well-known shape invariant 
Hamiltonians of the radial oscillator and the Darboux-Poschl- Teller potential to the 
shape invariant potentials of Odake-Sasaki. Dutta and Roy derived the two lowest 
members of the exceptional Laguerre polynomials by this method. The method is ex- 
panded to its full generality and many other ramifications, including the aspects of 
generalised Bochner problem and the bispectral property of the exceptional orthogonal 
polynomials, are discussed. 



1 Introduction 

Here we construct in an elementary way the four sets of infinitely many shape invariant 
Hamiltonians and the corresponding exceptional (Xt) polynomials [U |2]. The idea is quite 
simple. We start with a prepotential Wt(x\ A) and define a pair of factorised Hamiltonians 
(A) which are intertwined by the the Darboux-Crum transformations [31 H] in terms of 
Ai{\) and ^(A)t : 

H { +\\)=A\(\)M\), H<f\\) d ^MX)Al(\), £=1,2,..., (1.1) 



AeWU^W = H<r\X)M\), 4(X)H^(\) = U { e + \\)Ai(\), (1.2) 

. /.Ndef d dW e (x;X) t d dW e (x;X) 

A A A J = -, 3 j A?l A J = ~j j • v 1 - 6 ) 

ax ax ax ax 

Here A stands for the set of parameters of the theory. The prepotential We(x; A) is so chosen 
that / H^ r \X) is the well-known shape invariant [5] Hamiltonian of the radial oscillator (6j[7] 
potential or the Darboux-Poschl- Teller (DPT) potential [8| and \X) is the Hamiltonian 
of the recently derived shape invariant potentials of Odake-Sasaki [H |2j [9j [10] . The essential 
property of the prepotential to achieve the above goal is that e ±Wt ^ x ' x ^ is not square integrable. 
We know that T-L^\\) has a well defined groundstate. Thus these two functions e ±w t( x ' x ) 
cannot correspond to the groundstates of "H^^A) or f^~' ) (A). That is, the groundstates of 
?{^(A) are not annihilated by Af(X) or A\(X). This means that Uf\X) and U^ ] {X) arc 
exactly iso-spectral including the groundstates: 

n i ;\\)4^(x-,\) = s^(x)<f ) ^(x ] x) 

HjrHWftte A) = S^(X)4^(x; A) 



, 4t ) (A)=4; ) (A)>0, (1.4) 
1,2,..., 7i = 0,l,2,...,. (1.5) 



By construction "H^ (A) is shape invariant and exactly solvable. That is, the set of eigenval- 
ues {£e~ti{X)} and the corresponding eigenfunctions {(ft^Jix; A)} are exactly known. Through- 
out this paper we choose all the eigenfunctions to be real. They form a complete set of 
orthogonal functions: 

J <l>tm( x i X)(j) { ^{x] X)dx = h^ n (X)5 mn , h e , n (X) > 0. (1.6) 

Thanks to the intertwining relations (11.21) the eigenfunctions {<f>(^ (x; \)} of the partner 
Hamiltonian TL^~\X) are obtained from the well-known eigenfunctions {4*1^ ( x 'i X)} of the 
radial oscillator or the DPT Hamiltonian ^^(A) by the Darboux-Crum transformation in 
terms of Ag(X) 



>5; ) (x;A)=A £ (A)0(+ ) (x;A), £=1,2,..., n = 0,1,2,..., (1.7) 
and vice versa: 



iit ) (x;A) = 4^T^"'' 

4: } (A) 



(^ A ) = ^T77tC(^ A )' ^=1,2,..., n = 0,l,2,...,. (1.8) 
Of course {<p^(x; A)} is another complete set of orthogonal functions: 



4;l^x)4;n(x;x)dx = J mx)<^(x ] \)'M>)<&^^ 



= / 02(x;A)-^(A)%(A)0W(x;A) 

= 4t } W / A)d* = £^{\)h^{X)5 mn . (1.9) 

These imply the completeness of the exceptional orthogonal polynomials and the above 
relationship (jl.7p provides the formula relating the exceptional orthogonal polynomials to 
the classical orthogonal polynomials (i.e, the Laguerre or Jacobi polynomials) as shown in 
(2.1) and (2.3) of [10J. The orthogonality (jl.9p corresponds to the integration formulas 
derived in §7 of [10J. These will be demonstrated in detail in subsequent sections. 

The above requirements lead to the following general form of the prepotential Wi(x\ A): 

W t (x) A) = w (x; X + £S)+ log&Mx); A), (1.10) 

in which wo(x; A) is obtained by changing the sign of one term of the prepotential wo(x; A) 
corresponding to the radial oscillator or the DPT potential. Here 8 is the shift of the 
parameters. The change of the sign ensures the non-square integrability of e ±We{x < x l The 
additional term is the logarithm of the degree I eigenpolynomial, the Laguerre or Jacobi 
polynomial with twisted parameters or arguments, introduced by Odake-Sasaki [TJ [21 [10] . 

Before going to the details in the subsequent sections, let us make a few remarks about 
the background. This type of approach of deriving a new exactly solvable Hamiltonian from 
a known one in terms of Darboux-Crum [21 H] transformations has a long history and various 
aspects [HI [121 US]- The method we are concerned in this paper is, in its essence, based 
on an alternative factorisation of an exactly solvable Hamiltonian (plus a constant), e.g., 
the radial oscillator and the DPT potential. Some refer to those newly found Hamiltoni- 
ans as "conditionally exactly solvable". Junker and Roy [13] discussed an example of an 
alternative factorisation of the radial oscillator Hamiltonian by using the confluent hyper- 
geometric function iF±, which could encompass the results of the LI exceptional orthogonal 
polynomials [1] if the parameters and settings are properly chosen. After the introduction 
of the X\ Laguerre polynomials by Gomez-Ullate et al [H] and Quesne [151 02] an d the Xe 
(£ — 1, 2, . . .) polynomials by Odake-Sasaki [H [2], Roy and his collaborator [17] derived the 
X\ and X2 Laguerre polynomials of the LI type in this way. This will be mentioned in a 
later section. A recent report by Gomez-Ullate eta al [TS] has small overlap with the present 
work. So far as we are aware of, a "conditionally exactly solvable" treatment of the fully 
general (non-symmetric) DPT potential does not exist. Therefore, the present derivation of 
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the Jl and J2 exceptional Jacobi polynomials from the classical Jacobi polynomials based 
on the Darboux-Crum transformation is new. 

The exceptional orthogonal polynomials were originally introduced in [T3] by extending 
Bochner's theorem |19] for Sturm-Liouville problems. The characterisation of the exceptional 
orthogonal polynomials as polynomial solutions of Sturm-Liouville type equations under gen- 
eralised Bochner problems will be discussed in section HI We will show that the exceptional 
Laguerre and Jacobi polynomials have the bispectral property [20]. On top of the Sturm- 
Liouville type equation, the X^ polynomials satisfy A£ + 1 recursion relations (I4.15P and 
(I4.33p . which could be rewritten as an eigenvalue equation ( I4.16P of a semi- infinite matrix 
K which acts on the label n of the exceptional polynomials P n (x) and L n (x). Simple inter- 
pretation of the bispectral property is provided as the characteristic feature of the invariant 
polynomial subspaces of the Sturm-Liouville type operator. 

This paper is organised as follows. In section |2j the two types (LI and L2) of infinitely 
many exceptional Laguerre polynomials are derived by simple Darboux-Crum transforma- 
tions connecting them with the Laguerre polynomials. The two types (Jl and J2) of infinitely 
many exceptional Jacobi polynomials are derived in a similar way by using Darboux-Crum 
transformations connecting them with the Jacobi polynomials. The aspects of the gener- 
alised Bochner problems, in particular, the bispectral property will be discussed in section 
HJ Based on the so-called Darboux transformations of the orthogonal polynomials, that 
is the Christoffel and Geronimus transformations of the Jacobi and Laguerre polynomials 
[2T| [22], the bispectral property of the exceptional polynomials are derived elementarily. 
In other words, the Xi polynomials are shown to satisfy A£ + 1 term recurrence relations, 
which are the generalisation of the well known three term recurrence relations. The final 
section is for comments and discussions. It provides a simple proof of shape invariance of 
the Hamiltonians of the exceptional orthogonal polynomials. It is shown that they also 
possess the creation/annihilation operators inherited from the radial oscillator/DPT Hamil- 
tonian systems. It is stressed that the Hamiltonians of the radial oscillator/DPT potentials 
admit infinitely many non-singular factorisations, which induce the Darboux-Crum trans- 
formations. The extended 'three term recurrence' relations for the exceptional orthogonal 
polynomials are also discussed from the Darboux-Crum transformations point of view. 
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2 Exceptional Laguerre polynomials 



Here we will derive the LI and L2 exceptional Laguerre polynomials as well as the corre- 
sponding Hamiltonians, that is the potentials. 

2.1 Radial Oscillator 

Let us start with the radial oscillator with A = g > and 5 = 1: 

w (x; g) = f — \x 2 + glogx, < x < oo, (2.1) 

H%\g) = p 2 + x 2 + 9 -^± -1-2*. (2.2) 

x z 

It is trivial to verify the shape invariance [5]: 

n { R \g)=>H { +\g + l)+4, S 1 (g)=A. (2.3) 

Its eigenvalues and eigenf unctions are 

£n(g) = 4n, (2.4) 
n (x;s) = P„(r^)e^), 77 = r/(x) = x 2 , P n (x;<7) = pJT^x), (2.5) 

in which iffl (x) is the Laguerre polynomial satisfying the differential equation 

xd 2 x L^ (x) + {a + 1 - aOo^LW (x) + nL^ (x) = 0. (2.6) 

It should be stressed that the groundstate wavefunction </>o(x; g) = e w °( x '' g ^ = e~ x2 ^ 2 x 9 is 
square integrable and it provides the orthogonality measure of the Laguerre polynomials; 

e 2 Wo{x ; g)Pm{x 2. g)Pn{x 2. g)dx = hn ( g )S mn , h n (g) = f T(n + g + ±). (2.7) 
o 2n\ 2 

The radial oscillator Hamiltonian system is exactly solvable in the Heisenberg picture, too 
[23] . The exact annihilation/creation operators are obtained as the positive/negative energy 
parts of the Heisenberg operator solution (see, for example, (3.8) of [23J): 

a« = ((l,*f-^)/*. (2.8) 
The action of these operators are 

a { - ] (j) n (x] g) = -(n + g- -)</> n _i(x; g), a {+) <p n (x] g) = -(n + l)0 n+1 (x; g). (2.9) 
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2.2 LI and L2 exceptional Laguerre polynomials 

Here we derive the LI and L2 exceptional Laguerre polynomials. For each positive integer 
£ = 1,2,..., let us consider the pair of Hamiltonians (<?) and 7l\ \g) corresponding to 
the following prepotentials (rj = i](x) = f x 2 ) 

2 

LI: Wt(x;g) = y + (<? + £- 1) logx + log&(^); ^ g > 1/2, (2.10) 

£^)= f 4 9+ ^ f) (-^ (2.H) 

L2: ^(^^^-y-^ + ^logx + log^^^);^), g > -1/2, (2.12) 

^;S) = 4~ 9 ~ M) (^ (2-13) 

?4 +) G?) = ^(ff)^(ff), = ^(0)4(0), (2-14) 

^ = d _ MVpgl {215) 

By simple calculation using the differential equation for the Laguerre polynomial (|2.6p . the 
Hamiltonian (g) is shown to be equal to the radial oscillator with g — » g + £ — 1 for LI 
and with g — > g + £ + 1 for L2 up to an additive constant: 

LI: n^(g)=p 2 + X 2 + i9 + £ - 1 ^ 9 + i - 2) +2g + 6l-l, (2.16) 
= ?4 +) G? + * - 1) + 2(2^ + 4£ - 1), (2.17) 

L2: ft< + >(g) = P 2 + x 2 + ( ^ + i){9 + i+1) +2(g-£)- 1, (2.18) 

or 

= U { +\g + £ + !) + 2(2^ + 1). (2.19) 
The partner Hamiltonians are 



LI: =p 2 + a; 2 + v ^ ' "'^ 2 - ^ + 2{g - £) - 3 

\Umg) J V x + Uv;g) J' { } 

L2: Ut\g)=p 2 + x 2 + i9 + i){9 2 +i - 1) +2(g + 3£) + l 

x z 

n ( dMmg) \ ( q f | (g + t) \ d x Uv;g) \ (99U 
+ VUmg) A V x 7 fctos) 7' ( j 

Up to additive constants, the above Hamiltonians ( 12 . 20 j) and ( I2.2ip are equal to the Hamilto- 
nians of the LI and L2 exceptional orthogonal polynomials derived by Odake-Sasaki [TJ ITU] : 

n?~ s (9) =V + ( dw j X > 9) ) 2 + (2.22) 
V ax 1 ax 1 



(g+e)( g+ e-i) 



(/) = -— + (g + i) logx + log ^ , (2.23) 

LI: H < f\g)='H?- s (g) + 2{2g + U-1), (2.24) 
L2: H<f\g) = H?- s (g) + 2(2g + l). (2.25) 

The definition of £,e(r};g) for the LI and L2 Odake-Sasaki cases are the same as those given 
in (ETTT]) and fl2~T3l . 

Let us note that e Wi ^ x ' 9 ' is not square integrable at x = oo and e _H/<?(a:;9 ) is not square 
integrable at x = for the LI case, whereas for the L2 case e w ^ x;9 > is not square integrable 
at x = and e~ Wl ^ x '' 9 ^ is not square integrable at x — oo. In both cases the prepotential 
Wi(x; g) ( 12.1 Oft and ( 12. 12ft are regular in the interval < x < oo. For this, it is enough to 
show that £i(j)(x); g) does not have a zero in < x < oo. In fact we have, 

Ll: ^y,9) = j: {S+ ^ e k _~f'- t ^>0, (2.26) 



k=0 

L2: (-l)«^W; S ) = ^|±i!ii^>0, (2.27) 

as shown in (2.39) of [9J. Thus we find that e ±We ^ x ' 9 ^ cannot be the groundstates of the 
Hamiltonians (<?)■ However, we know quite well that "Hi (g), being the radial oscillator 
Hamiltonian, has a well-defined groundstate. This means that the partner Hamiltonian 
T-L[~\g), thus the Hamiltonian of the Ll and L2 exceptional Laguerre polynomials, are 
exactly iso-spectral to the radial oscillator Hamiltonian 7/^(p), which have the following 
eigenvalues and the corresponding eigenfunctions: 

Ll: 4t ) (o) = 4n + 2(2o + 4£-l), 0«(x; g) = e'^x^L^ } (n), (2.28) 
L2: S^\g) = An + 2(2o + l), <${x;g) = e-^x 9+£+1 L ( ^\rj). (2.29) 

The intertwining relation (jl.2p implies the simple expressions of the eigenfunctions of the 
partner Hamiltonians %^ in terms of A^g): 

Ll: 4; ) (o) = 4n + 2(2o + 4£-l), £° n s (o) = 4n, (2.30) 



L2: e< c ;J(3) = tn + 2(2g+l), £ °f( 9 ) = in, (2.32) 



^fo^W^^fo)). (2.33) 



These are to be compared with the explicit expressions of the exceptional Laguerre polyno- 
mials Pe, n {T]',g) derived in [TO]: 



<Pe, n (x; g) = il>i(x\ g)Pt, n (rj; g), M x ; g) = e c , X \ , 9) = \ L i a ^k v } ' L1 

Uv,g) { L\ 9 2 \t]) : L2 



^»fa*) = { (n + . + I^^ + I^. + l)^^^) • (0-S2.1) 

+ ^(r/;^)^Li 9+£+|) (r ? )) : L2, 

The final expressions of the eigenfunctions (I2.3ip and fl 2 . 3 3 j) are the same as (2.1) of [10] up 
to a multiplicative constant. Use is made of the identities of the Laguerre polynomials (E.ll) 
and (E.12) of [ID] . Thus we have derived the Hamiltonians as well as the eigenfunctions, that 
is, the LI and L2 exceptional Laguerre polynomials and the weight functions, from those of 
the radial oscillator by the Darboux-Crum transformations. 

3 Exceptional Jacobi polynomials 

Here we will derive the Jl and J2 exceptional Jacobi polynomials as well as the corresponding 
Hamiltonians, that is the potentials. 

3.1 Trigonometric DPT potential 

The trigonometric DPT [S] potential has two parameters A = (g, h), g > 0, h > and 



1,1 , 



wq(x; g, h) = f g log sin x + h log cos x, < x < — , (3.1) 

Kdpt(9, h)=p + —r- 2 — + 2 (g + h) . 3.2 

sin £ cos^x 
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It is trivial to verify the shape invariance [5]: 

H ( dIt(9, h) = n^ PT (9 + l,h + l) + 4(g + h+ 1), E x {jg, h) = 4(g + h + 1). (3.3) 

Its eigenvalues and eigenf unctions are: 

£ n (g,h) =4n{n + g + h), (3.4) 

<t> n {x ] g,h) = P n {r ]] g,h)z Wo{x '' 9 ' h \ V = v(%) = cos 2x, P n (x; g, h) = P^' h ~ l \x), (3.5) 

in which Pn (x) is the Jacobi polynomial satisfying the second order differential equation 

(1 - x 2 )d 2 x Pi a ^ (x) + (P-a-(a + P + 2)x) d x P^ ] (x) +n(n + a + + 1)P^ (x) = 0. (3.6) 

The groundstate wavefunction (f>o(x;g,h) = e Wo ^ x ' 9 ' h > = (sina;) 9 (cosa;) /l is square integrable 
and it provides the orthogonality measure of the Jacobi polynomials: 

»tt/2 

e 2m ^ sM P m {r,(x);g,h)P n {ri(x);g,h)dx = h n (g,h)S mn , (3.7) 

vy ' 7 2n!(2n + # + /i)r(n + # + /i) v ; 

The trigonometric DPT potential is also exactly solvable in the Heisenberg picture [23]. The 
annihilation and creation operators are (see, for example, (3.28) of |23j): 

a'W/2 = a [±) 2VU' = ± sin 2x-$~ + cos 2x VU' + 01 ~ @ (3.9) 

dx y/W±l 

in which 

U"^U DPT +(g + hf, a^g- 1 -, P^h- 1 -. (3.10) 

When applied to the eigenvector <p n (13. 5 p as £ n (g, h) + (g + h) 2 = (2n + g + /i) 2 , we obtain 
(see, for example, (3.29) and (3.30) of [23J): 

a 1 ; 0„(x; g>, h)/2 = — sin2x h (2n + # + h) cos2x0„(x; g, a) 

2n + a + p 
A{n + a){n + (3) 

= — 7 — 7~1T~ ^n-iix; g, h), 3.11 
2n + a + p 

+ a — 1 1 a 1 o ^« x; f ' ^ 

2n + a + p + 2 

4(n + l)(n + a + /3 + l) 

= 2n + a + /? + 2 ( 3 - 12 ) 



3.2 Jl and J 2 exceptional Jacobi polynomials 

Here we derive the Jl and J 2 exceptional Jacobi polynomials. As explained in [9j [10], the 
exceptional Jl and J2 orthogonal polynomials are 'mirror images' of each other, reflecting 
the parity property Pi a '^ of the Jacobi polynomial. Here we present 

both cases in parallel so that the structure of these polynomials can be better understood 
by comparison. 

For each positive integer £ = 1,2,..., let us consider the pair of Hamiltonians (g, h) 
and T-Li~\g, h) corresponding to the following prepotential (rj = r](x) = f cos2x) 

Jl: Wg(x; g, h) = f (g + £ — 1) log sin x — (h + £) log cos x + log ^(77; g, h), (3.13) 

U^gM^P^^-^^), g>h>0, (3.14) 
J2: We(x; g, h) = — (g + £) logsinx + (h + £ — 1) logcosx + log<^(?7; g, h), (3.15) 

Umg,h) = Pi~ g ~ i ~* Mt ~*\'l)> h>g>0, (3.16) 
Hf\g,h) ^A\(g,h)Ae(g,h), H { f\g,h) = f A e (g , h) A\(g , h) , (3.17) 

A( , v def d dW £ (x;g,h) t d dW e (x;g,h) 

Mg> h ) = ~j ~j > A\{g,h) = ~j • ( 3 - 18 ) 

By simple calculation using the differential equation for the Jacobi polynomial (13. 6p . we 
obtain the trigonometric DPT potential for 1-L^\g, h) up to an additive constant: 

Jl: nf\ S M=P> + {9 + l - 1){ ? + e - 2) + ih + e)i ^ e+1) -(2i + S -h-l)\ 

sin x cos^ x 

(3.19) 

= ?4 + pt(# + £-l,h + £ + l) + (2g + 4£-l)(2h+ 1), (3.20) 

J2: ««( 9 , h) = P 2 + fe + W+D + (ft + <-D(* + <-2) -(« + *-,- 1)=, 

sin x cos z x 

(3.21) 

= ^D } p T (g + £+l,h + £-l) + (2h + 4£-l)(2g + l) (3.22) 
The partner Hamiltonians are 

Jl: H^,h)=p' + {3 + li{3 , + t ~ 1) + lh + li{ \ +l ~ 1) 

sin x cos z x 

+ 2 ( d : ^9,h) \ ( 2{{g + £ _ 1)cotx + {h + e)t&nx) + d Mv,g,h) 



&(v;g,h) J V £t(v;g,h) 

~{2£ + g-h-l) 2 + 8£{£ + g-h-l), (3.23) 

j 2: ««-»(,,/,) =p . + to + Qto + <-i) + (* + Qfr + <-i) 

sin a; cos" 1 x 
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+ 2 [ d : UV;9 : h : ) ) (-2 {{g + £) cot x + (h + i- 1) tan*) + 9 ^ 9 > k) 



-{2£ + h-g-lf + U(£ + h-g-l). (3.24) 

Up to additive constants, the above Hamiltonian f!3.23j) and f)3.24p are equal to the Hamil- 
tonian of the Jl and J2 exceptional Jacobi polynomials derived by Odake-Sasaki [Tl [TO]: 



H?- s (g,h)-P 2 +( dMX ; 9,h) ) + * We< ?[ 9 > h \ (3.25) 
\ dx J dx 2 

we(x; g, h) = f (g + £) log sin x + (h + I) log cos x + log ^—h - — — (3.26) 

Jl: HPfah) =U?- S (g,h) + (2(7 + 4£ - 1) (2/i + l), (3.27) 
J2: V!f\g,h) = n^ s (g,h) + (2h + 4£-l)(2g + l). (3.28) 

The definition of £,e{rj; g) for the Jl and J2 Odake-Sasaki cases are the same as those given 
in (13.141) and (I3.16p . Again let us note that e Wt ^ x > 9,h ' is not square integrable at x = n/2 and 
e -w e (x;g,h) j g nQ ^ S q Uare integrable at x = for the Jl case, whereas for the J2 case e We ^ x ' 9 ' h ^ 
is not square integrable at x = and e~ We ^ X]9 ' h,> is not square integrable at x = ir/2. But the 
prepotentials Wg(x; g, h) ( 13 . 1 3[) and (I3.15P are regular in the interval < x < ir/2. For this, 
it is enough to show that £,e{vi x )'i g, h) does not have a zero in < x < it/2. In fact we have, 

Jl: (-ifefaM; 9, h) = ^ t " " Wrf T^Pl = 1> " > °- < 3 ' 29 > 

« ! k\(h + £ - k + i) fc 

J2: (- lfe („ W ; g , h) = ^ E ^ *,+ 'jf J^'." 1)t »)" > 0, (3.30) 



fc=0 



as shown in (2.40) of [9]. Thus we find that e ±We ^ X]9 ' h ^ cannot be the groundstates of the 
Hamiltonians "H^ h). However, we know well that 7i^\g,h), being the trigonometric 
DPT Hamiltonian, has a well-defined groundstate. This means that the partner Hamiltonian 
f-L^ig, h), thus the Hamiltonian of the Jl and J2 exceptional Jacobi polynomials, are exactly 
iso-spectral to the trigonometric DPT Hamiltonians T-i\ <(g,h), which have the following 
eigenvalues and the corresponding eigenfunctions: 

Jl: £$(g,h) =An(n + g + h + 2£) + (2g + A£-l)(2h+l) J (3.31) 
4^(x;g,h) = (sinxy +e -\co S x) h+e+1 Pi 9+e - hh+e+l2) ( V ), (3.32) 
J2: £i + n \g,h) =4n(n + g + h + 2£) + (2h + A£-l)(2g + l) } (3.33) 
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(3.34) 



The intertwining relation (11.2j) implies the simple expressions of the eigenfunctions of the 
partner Hamiltonians T-i\ in terms of Ae(g): 

Jl: Si~\g,h) = 4n(n + g + h + 2£) + (2g + A£-l)(2h + l), (3.35) 
£°; s (#,/i) = 4n(n + # + /i + 2£), (3.36) 



cf>°^(x-,g,h) = ( f>\J(x;g,h) 



(-)< 



-2 



d dWe(x;g,h) 
doc dec 

x (sinx)^ £ - 1 (cosa;) h+£+1 Pi !;+£ " f ' ,l+ ' +|) (r ? ) 



(v) 

(3.37) 
(3.38) 
(3.39) 



J2: ^(p, h) = An(n + g + h + 21) + (2h + U - 1)(1 + 2#), 
^°; s (^,/i) =4n(n + (? + /i + 2£), 

x (sinx) 5+£+1 (cosa;) /l+£ - 1 Pi 9+£+i ' l+ '" |) (77) 

(sinx)^(cosx) w / 1 (o^+i^-D 
= -2 ^7 [{9 + :;Kmg + l,h + l)Pn {V) 



-(1 - vMv,9,h)d v P!f +i+l >' h+e ^ \ V ) ) . 



(3.40) 

These are to be compared with the explicit expressions of the exceptional Jacobi polynomials 
Pe,n(v'i 9 ') h) derived in [TO]: 

, / , N , / , \ n / ,\ , / , \ def (sinx) 9+e (cos x) h+t 
(pe,n{x; g, h) = ip e (x; g, h)P e Jr]; g, h), ipi{x; g,h) = - 



PeAv,9,h) = < 



&{v;g,h) 

Uu a h) * I > * > : Jl 

{ P £ 2 2 (r/), h> g > : J2 

(n + /i + J)" 1 ((A + <? + 1, + l)P? + '~ f , " +£+ ^ ) (7 7 ) 

+ (1 + vMm9,h)d r ,Pi 9+i ^' h+i+ ^(r 1 )) 

(n + g + I)- 1 ((g + g + i,h+ i)P^' h+ ^\ v ) 



: Jl 



(0-S2.3) 



The final expressions (I3.37|) and ( I3.40p are the same as (2.3) of [10] up to a multiplicative 
constant. Use is made of the identity of the Jacobi polynomials (E.22) of [10]. Thus we have 
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derived the Hamiltonians as well as the eigenfunctions of the Jl and J2 exceptional Jacobi 
polynomials from those of the trigonometric DPT by the Darboux-Crum transformations. 

The exceptional Laguerre and Jacobi polynomials satisfy a second order linear differential 
equation in the entire complex rj plane: 



nf- s {X)P e , n {m A) = S e!n {X)P tin { m A), £ An (A) = E n {\ + id). 



(3.41) 



For later use we give the explicit form of the second order Fuchsian differential operator Hf~ s 
which was given in (3.5) of [10J: 



^°" S (A) = -4(02(77)^ + ( Cl ( V ,X + £5) -2c 2 (r))d v logU m X))j- 



drj 

w) 6w;A) 4 

in which ci, c 2 , d\, d 2 and £e are given by 
def / # + 5 — 77 :L 



(3.42) 



ci(v, A) 



h-g-(g + h + l)r] :J 



def J 77 : L 

1 - 77 2 : J 



^(A) 



def 



9 + | 



LI 

L2,J2 
Jl 



d 2 (v) = \ -'I 



LI 
L2 

J1/J2 



def f =F4£ 



: L1/L2 



A£(£±gTh-l) : J1/J2 



(3.43) 
(3.44) 
(3.45) 



4 Generalised Bochner problem: bispectral property 

The exceptional Laguerre polynomials LI and L2 as well as the exceptional Jacobi polynomi- 
als Jl and J2 belong to complete orthogonal families of functions (the completeness follows 
from the well known properties of the Darboux process; indeed, the Darboux transformation 
which does not generate new eigenstates preserves the completeness of transformed system 
of eigenfunctions as solutions of the self-adjoint Schrodinger equation). These polynomials, 
however, do not belong to the ordinary families of orthogonal polynomials because polyno- 
mials of the first i — 1 degrees are absent in these systems [14J. Hence these polynomials 
do not satisfy 3-term recurrence relation which is a characteristic property of nondegener- 
ate orthogonal polynomials (see, e.g. |24j). Nevertheless, as we will show, the exceptional 
polynomials Jl, J2, LI and L2 do satisfy A£ + 1-term recurrence relations, i.e. they can be 
considered as eigenvectors of a semi-infinite matrix K having 4£+l diagonals. In this sense 
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the considered exceptional polynomials possess a very important bispectral property [20J: 
they are simultaneously eigenf unctions of a Sturm- Liouville operator and a matrix K. 



4.1 Bispectrality of the exceptional Jacobi polynomials 

We consider exceptional orthogonal polynomials of the J 1- type (the type J2 can be considered 
in the same manner). For simplicity of presentation we slightly change the previous notation, 
denoting the exceptional Jl polynomials as P n (x) and introducing standard parameters 



a d =? g + I - §, b = h + £ + i. We also use x for the sinusoidal coordinate 77 and denote 
tt(x) = £,e(x; g, h) = P^ a ' b \x) which is a polynomial of degree £. This polynomial will play a 
crucial role in the following. Then formula (0-S2.3) for the Jl polynomials can be presented 
in the form 

P n (x) = tt(x) ((1 + x)P% a > b \x) + bP^ b \x)) - tt / (x)(1 + x)P^ h \x), (4.1) 

which is equal to P^ n (x; g, h) in (0-S2.3) up to a multiplicative constant. Recall that the 
Jacobi polynomials Pi a ' b ^ (x) satisfy the three-term recurrence relation 

AnPt?(x) + B n P^ b \x) + C n P^{x) = xP^ b \x), (4.2) 



with 



2(n + l)(n + a + b + 1) b 2 - a 

An — ~T~~ : rr~7Z : TT , D 



(2n + a + b+l){2n + a + b + 2Y " (2n + a + b)(2n + a + b + 2) 
2{n + a){n + b) 



(2n + a + b)(2n + a + b + 1) ' 
They are orthogonal on the interval [—1, 1] 



P^ b \x)P^ b \x){l - x) a {l + xfalx = h n S nm , (4.3) 



where 

, dc f 2 a+b+1 T(n + a + l)T(n + b + l) 

K 



n\(2n + a + b+ l)T(n + a + b + 1) 
is the normalization constant. The exceptional Jl polynomials are orthogonal on the same 

interval 

(P n , P m ) = J P n {x)P m {x)w{x)alx = h n 5 nm , (4.4) 
with the weight function 

w{x) = ^— 4.5 
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and some nonzero normalization coefficients h n (in fact, these coefficients can easily be 
connected with the coefficients h n , however we need not their explicit expressions here). 
Using elementary properties of the Jacobi polynomials [21], or (2.24) of [ID] , we also have 

P n (x) = (b + n)7r(x)Pi a+1 > b - 1) (x)-ir'(x)(l + x)Pi a ' b) (x). (4.6) 

Moreover, there are classical formulas [2T] 

(2n + a + b + 1)(1 + x)P^ b \x) = 2{n + l)P^\~ l \x) + 2(n + b)P^ h - x \x) : (4.7) 

and 

(2n + a + b + l)Pi a ' fe) (^) = (n + a + 6 + l^+^x) - {n + b) P^X^ (x) . (4.8) 

These formulas have a simple interpretation in terms of so-called Darboux transformations 
of the orthogonal polynomials [22]. Namely, the formula (14. 7p describes the Christoffel 
transformation of the Jacobi polynomials P^ ,b ~ l \x) — > Pn a ' b \x) while the formula (14. 8[) 
describes the Geronimus transformation Pn (X). 
Combining formulas (" 14. T[) and (14. 8p we obtain 

il + x)P^ b \x) = o^P^%)+^Pi^ 1 \x)+ 7n P^%), (4.9) 

def 2(n + l)(n + a + b+ 1) de f 2(a + &)(« + &) 

{2n + a + b + l)(2n + a + 6 + 2) ' ^ n (2n + a + 6)(2n + a + 6 + 2) ' 
def 2(ti + 6)(ti + 6-1) 
7 " ~ {2n + a + b){2n + a + b+ 1) ' 

Using the formula (14. 9 p with the repeated use of the three term recurrence relation ( 14. 2 p for 
the factors tt(x) and tt'(x) in (14.61) . we conclude that 

n+l 

Pn(x) = ^sP^-'Kx), (4.10) 

s=n— I 

with some real coefficients 

The formula (14.101) establishes a "local" property of the exceptional polynomials P n (x) 
with respect to the basis Pn a+1 ' b ~ 1 \x), n = 0,1,2,.... Recall that deg(P n (x)) = n + i 
and hence degrees of polynomials in lhs and rhs of (I4.10p coincide. There is a reciprocal 
"local" relation which expresses polynomials tt 2 (x) Pn a+1 ' b X \x) in terms of a finite linear 
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combination of the exceptional polynomials P n {x). In order to get this relation, we can 
expand the product of polynomials tt 2 (x) Pn +1 ' b (x) in terms of the basis P n (x): 

oo 

7T 2 (X) Pi^^ix) = J2VnsP S (x). (4.11) 

s=0 

This sum in general can be infinite, because the polynomials P n (x) form a basis in a Hilbert 
space with the scalar product given by the formula (14.41) but due to existence of a "gap" in 
degrees of polynomials P n (x) for a generic polynomial Q(x) we should have an expansion 

oo 

Q{X) =J2(nsPs{x) 
s=0 

with infinitely many coefficients ( ns . However, for some special choices of the polynomial 
Q(x) this expansion can contain only a finite number of terms. 

In order to find the coefficients rj ns , let us multiply both sides of (14. lip by w(x)P s (x) 
(w(x) is given by (14.51) ) and integrate over the interval [—1, 1]. Due to the orthogonality 
relation for polynomials P n (x), in the rhs after integration we obtain the term h s rj ns with 
the nonzero coefficient h s in (14.41) . On the other hand, in lhs we have the integral 

• i 

(1 - x) a+1 (l + x) b - 1 Pj i a+1 ' b ~ 1 \x)P s (x)dx. (4.12) 



-i 

Substituting expression (I4.10p into (14. 12j) and using the orthogonality property of the Jacobi 
polynomials Pn +1,b ~ 1 \x) we have the relation 

h s 7] ns = h n £ sn . (4.13) 

Hence only H + 1 coefficients r] ns are nonzero. We thus have the expansion 

n+e 

-k\x) P^ b - l) {x) = Vn,P s (x), (4.14) 

s=n—i 

where the coefficients r] ns are connected with £ ns by the "mirror" relation (14.131) . Formulas 
(I4.10p and (14.141) can be considered as a generalization for generic I of corresponding formulas 
obtained for i = 1 in |14j . 

Eliminating the Jacobi polynomials Pn +1 ' b ^(x) from these formulas, we arrive at the 
recurrence relation 

n+2l 

7l 2 (x)P n (x)= K nsPn(x), P s (x) = 0, if 8 < 0, (4.15) 

s=n~2i 
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with some real coefficients K ns . This recurrence relation belongs to the class of 4£+l-diagonal 
relations. This means that in the operator form the relation (14. 15j) can be presented as 

KP = tt 2 (x)P, (4.16) 

— * 

where P is an infinite dimensional vector with components P n (x), n = 0, 1, 2, . . . , 

P = {P (x),P 1 (x),P 2 (x),...}, 

and K is a matrix with entries K ns . This matrix has no more than 4£+l nonzero diagonals. 
Corresponding polynomials P n (x) satisfy 4£+l-term recurrence relation (|4.15p . The ordinary 
orthogonal polynomials satisfy 3-term recurrence relation. Hence we have polynomials P n {x) 
satisfying more general recurrence relation. Recurrence relations of such a type were studied 
e.g. by Duran and Van Assche [25] who showed that such polynomials should satisfy a 
matrix orthogonality relation. Thus the exceptional Jacobi polynomials belong to the class of 
polynomials satisfying higher-order recurrence relations. Note however, that in the approach 
of [25] it is required that polynomials P n (x) have exactly degree n = 0,1,2, ... . In our case 
the polynomial P n {x) has degree n + i, n = 0, 1, 2, . . . . This means that the methods of [25] 
should be modified if applied to the case of the exceptional polynomials. 

Formulas ( I4.10p and ( I4.14p admit an interesting algebraic interpretation. Introduce semi- 
infinite matrices H and H by their entries £ ns and r] ns . Then we have 

tt 2 (x)P(x) = HEP(x), (4.17) 

where 

P(x) = {P^ +1 ' b - 1 \x),Pi a+1 ' b ' 1 \x),Pt +1 ' b - 1 \x),...h 

and 

n 2 (x)P = EHP, K = EH. (4.18) 
From the relation ( I4.17P it follows that 

tt 2 (J) = HE, (4.19) 

where J is a Jacobi (3-diagonal) matrix corresponding to the Jacobi polynomials P^ +l,h (x), 
i.e. 

xP{x) = JP{x). 
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Thus the matrices H and S appear under factorization of the 4£ + 1-diagonal matrix ir 2 ( J). 
The exceptional polynomials P n (x) satisfy recurrence relation (14. 18j) which appear after refac- 
torization (permutation) of the factors H and S. Such permutation of the matrix factors 
is known as the Darboux transformation of the matrix 7r 2 (J). Similar Darboux transfor- 
mations were already studied in [20] where refactorization of the quadratic polynomials in 
J corresponding to the Jacobi and Laguerre polynomials was considered. Clearly, under 
such refactorization one obtains new polynomials satisfying a 5-term recurrence relation. 
These polynomials are not classical orthogonal polynomials. Nevertheless, the authors of 
[20] showed that these new polynomials are eigenfunctions of a linear fourth-order differen- 
tial operator. In our case we have almost the same construction as in (20] but the resulting 
polynomials are exceptional (i.e. some degrees are absent) and satisfy a second-order differ- 
ential equation. 

It is also interesting to note that the weight function w(x) given by fj4.5j) . differs from 
the "classical" weight function (1 — x) a+1 (l + x) b_1 for the Jacobi polynomials by the factor 
7r~ 2 (x). For the ordinary orthogonal polynomials it is known that a rational modification 
of the weight function w(x) — > w(x)U(x)~ l with some polynomial U(x) of degree M is 
equivalent to the application of M Geronimus transforms [22J. If P n (x) are orthogonal 
polynomials corresponding to the weight function w(x) then orthogonal polynomials P n (x) 
corresponding to the weight function w(x)U(x)~ 1 are given by the linear combination [22J 

P n (x) = A^P n (x) + A^P n ^(x) + ■■■ + 4 m )P„_ m (x), (4.20) 

with some coefficients A„ ■ , % = 0, 1, 2, . . . M. The formula (I4.20p describes a special case of 
the M-time Geronimus transform (general case of the Geronimus transform allows an adding 
of M concentrated masses located on the roots of the polynomial U(x) |22j). 

If we now compare formulas (14.101) and (14.201) we see that the exceptional Jl Jacobi 
polynomials P n (x) are connected with the ordinary Jacobi polynomials Pn a+1 ' b (x) by a 
transformation which resembles the ordinary M-time Geronimus transform, where M = 21. 
However, this transformation can be considered as a degeneration of the ordinary Geronimus 
transform in a sense that the first I polynomials P%{x\ i — 0,1, ...,£ — 1 become zero. 

The case of the J2 exceptional Jacobi polynomials can be considered in the same manner 
leading to a recurrence relation of the same type (14.151) . 
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4.2 Bispectrality of the exceptional Laguerre polynomials 

Consider the case of the exceptional LI Laguerre polynomials. Here we list some well known 
formulas for the ordinary Laguerre polynomials (see, e.g. [2Tj): 

Recurrence relation: — (n + l)L^\i{x) + (2n + a + l)L^\x) — (n + a)L[ l a 2 1 (s) 

= xL^(x), (4.21) 
Differentiation formula: L^ ] \x) = -L^\x), (4.22) 

Geronimus transformation: (x) = L^ +1) (x) - (x) . (4.23) 

Following (12. lip , we have tt(x) = £e(x; g) = L^\—x), where a = f g + i — 3/2. We denote 
the LI exceptional Laguerre polynomials cis L n (x ), which are equal to Pe }n (x; g) given in 
(0-S2.1) up to an overall sign 

L n (x) = n(x)L^\x) - (tt(x) + it' (x)) L^(x) . (4.24) 

Here the Geronimus transformation (14.231) is used to rewrite £i(x; g + 1) in terms of 7r(x) = 
£t(x; g). These polynomials start from degree £, i.e. deg(L n (x)) = n + £, n = 0, 1, 2, . . . and 
there are no other members of this orthogonal polynomial family. The orthogonality relation 
reads 

POO 

(L n ,L m )= L n (x)L m (x)w{x)dx = h n S nm , (4.25) 
Jo 



where 



/£>Gi | 1 g X 



*(,) = -^ w . (4.26) 

Now we derive recurrence relations for the exceptional Laguerre polynomials. First of 
all, using (14.221) and ( 14.231) we rewrite the formula (14.241) in the form 

L n (x) = -v'{x)L^l\x) - {tt'{x) + n(x))Li a+1 \x). (4.27) 

The polynomial n(x) has degree I. Hence, by the recurrence relation (14.21 j) we have 

n+e 

Ln(x)= Zn S L^ +1 \x), (4.28) 

s=n—i 

with some real coefficients £ ns . This is a formula expressing the polynomials L n ( 
linear combination of the Laguerre polynomials L^ +1 \x). In order to obtain the reciprocal 
formula let us consider the expression 

oo 

n 2 (x)Li a+1 \x) = J2 Vnk L k (x), (4.29) 

fc=0 
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with some real coefficients f] n k- Multiply both sides of (I4.29P by L s (x) and integrate with 
the weight function (I4.26[) over the interval (0, oo). We then get 

POD 

h s Vns= / x a+1 e- x L^ +1 \x)L s (x)dx. (4.30) 
Jo 

Substitute the expression (I4.28[) for L s (x) in terms of L^ +l \x) into rhs of (I4.30p . Using the 
orthogonality relation for the Laguerre polynomials, we have 

h s 7] ns = h n £ sn . (4.31) 

We already have that £ ns — if \n — s\ > £. Hence r] ns = if |n — s\ > I and we have the 
expansion with 21 + 1 terms 

n+e 

tt 2 (x)L^ +1 \x) = VnkL k (x), (4.32) 

k=n— i 

where the coefficients r] ns are connected with £ ns by (14.3 ip . 

Similarly to the case of the Jl exceptional Jacobi polynomials we obtain a A£ + 1-term 
recurrence relation 

n 2 {x)L n {x)= K ns L s {x), L s {x)=0, if s < 0, (4.33) 

s=n-2l 

The case of the L2 exceptional Laguerre polynomials can be considered in a similar manner 
leading to a recurrence relation of the same type (14.331) . 

One can formulate a natural conjecture that all the exceptional orthogonal polynomials 
(i.e. polynomials satisfying a linear second-order Sturm-Liouville equation and orthogonal to 
one another, see the detailed description of this problem in [TJ]) satisfy a recurrence relation 
of the type (14.151) or (I4.33P with an appropriately defined polynomial ir(x). 

It is worthwhile to stress that the function ir(x) 2 plays the important role of the eigenvalue 
of the recurrence relation (I4.15P or (14.331) . In the ordinary three term recurrence relation 
(14. 2p or (14.211) . the same role is played by the (sinusoidal) coordinate x itself. 

4.3 Invariant polynomial subspaces 

Here we will show that the bispectral property of the exceptional Laguerre and Jacobi polyno- 
mials (14.151) and (I4.33P can be easily understood as the characteristic feature of the invariant 
polynomial subspaces of the second order Fuchsian differential operator "H°" s (13.411) . Because 
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of the denominators containing tt(x) = £i(x; A), the operator does not map a generic 

polynomial in x into a polynomial. In other words 

H?" S (A)V n % V n , V n = Spanfl, x, . . . , x n ], n = 0, 1, 2, . . . , . (4.34) 

However, it is easy to see that 

n?- s (\)7i(x) 2 V n C V„ +2 £, n = 0, 1, 2, . . . , . (4.35) 

In other words, 7r(x) 2 V n belongs to the degree n + 2£ invariant polynomial subspace of the 
operator T-Lf~ s , see §6 of [ID] . Therefore, 7r(a;) 2 P n (x) and 7r(a;) 2 L n (a;) also belong to the degree 
n + 3£ invariant polynomial subspace and they can be expressed as a linear combination of 

{P n (x)} ({L n (x)}y. 

n+2£ n+2l 

7r{x) 2 P n {x) = KnsP s (x), 7c{x) 2 L n {x) = ^A(s), n = 0, 1, . . . , . (4.36) 

s=0 s=0 

Next let us consider the inner products (I4.4p and (14.251) for s < n — 2£\ 

(rr(x) 2 P n , P s ) = (P n , n(x) 2 P s ), (rr(x) 2 L n , L s ) = (L n , n(x) 2 L s ). 

Since 

s+2t s+2i 

tt{x) 2 P s {x) = ^ K sjPj, rc(x) 2 L s (x) = ^ K sjLj, 

3=0 3=0 

we find 

(n(x) 2 P n ,P s ) = 0, (ir{x) 2 L n , L 8 ) = 0, K ns = 0, s<n-2£. (4.37) 
This simply means the desired results 

n+2l n+2£ 
7l(x) 2 P n (x)= K nsPs(x), 7T(x) 2 L n (x)= K nsL s (x), 71 = 0,1,...,. (4.38) 

s=n-2£ s=n-2i 

This is essentially the same argument for demonstrating the three term recurrence relation. 

5 Comments and discussion 

A few remarks are in order. Dutta and Roy [IT] derived the first two members of the LI 
exceptional Laguerre polynomials in a way shown in section |2j In our language, they used a 
prepotential 

dcf ^ 2 

Wt{x;g) = — + (g+ 1) logx + logu £ (x; g), (5.1) 
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^(x;g) d ^ f 4 9+|) (-x 2 ), £=1,2. (5.2) 

But erroneously they insist that the partner Hamiltonian T-L\ (g) is not shape invariant on 
account of the fact that its partner is the radial oscillator Hamiltonian. 

Let us emphasise that the shape invariance is the intrinsic property of the Hamiltonian, or 
the potential, determined by the unique factorisation in terms of the groundstate wavefunc- 
tion. Allowing for singularities, a generic quantum mechanical Hamiltonian H has infinitely 
many different factorisations. For example, if 

n = p 2 + V{x), H<t> n (x) = £ n n (x), (5.3) 

then it is trivial to show the following factorisation: 

f d d x (j) n (x)\ ( d d x cj) n (x)\ ^ ^ ^ 

V dx 4> n {x) ) \dx 4> n (x) J n ' _ " 

The right hand side as a whole is non-singular, but each factor is singular, since (f) n {x) has 
n zeros. Usually the groundstate n = provides the unique non-singular factorisation. 

The results of the present paper assert that the radial oscillator (DPT) potential ad- 
mits two families, LI and L2 (Jl and J2), of infinitely many non-singular factorisations. 
Perhaps it would be worthwhile to reflect upon the significance of the infinitely many non- 
singular factorisations in historical perspective. More than three decades ago, Miller [26] 
embarked on the program to classify and exhaust factorisations of shape invariant Hamilto- 
nians (H = —d 2 /dx 2 + V(x)). Although he did not employ the word 'shape invariance', the 
essential ingredients were there. He considers the cases in which the Hamiltonian depends 
on one parameter only, say, g. Based on an assumption, rephrased in our notation, that the 
derivative of the prepotential has a finite power dependence on g 

= ST a 3 (xW, N, M G Z + , (5.5) 

he came to the conclusion that the allowed factorisation types were the same as those listed 
by Infeld-Hull [6]. The present cases of the radial oscillator, for example the L2 f )2.12p 

dW t (x]g) 9 + ? d x &(ri; g) 
la: — —x i — — {O.o 

dx x ^{rj; g) 

are not covered by his assumption. 
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It should be remarked that the present method provides an alternative proof of shape 
invariance of the Odake-Sasaki Hamiltonians [TJ |2j [91 [10]. Here we will show the shape 
invariance of the Hamiltonians of the exceptional Laguerre polynomials. Starting from the 
prepotential Wi(x;g), (I2.10p for the LI and (12.121) for the L2, we obtain the Hamiltonian 
fl2TT6|) for the LI and fl2TT8|) for the L2 and the partner Hamiltonian H\~\g), fl2~20j) 
for the LI and (I2.2ip for the L2. Then we can construct another partner Hamiltonian of 
'H < i r \g) by using the factorisation in terms of the groundstate wavefunction e~ x2 / 2 x 9+e ~ 1 
for the LI and e - x2 / 2 x 9+e+l for the L2. By the shape invariance of the radial oscillator 
Hamiltonian, the result is simply 7i^\g + 1) + 4. Thus it is also obtained by the shifted 
prepotential We(x; g + 1), which produces the partner Hamiltonian h[ (g + 1): 

W e (x;g + 1) 



uf\g + i] 



n^\g + i) 



groundstate 
wavefunction 



groundstate 
wavefunction 



n { e +) (g) 



W e (x; g) 



n { e -\g) 



Laguerre 
polynomial 



LI, L2 exceptional 
Laguerre polynomial 



Two ways of proving shape invariance. 

The direct proof of shape invariance of ^^(g) by using the groundstate wavefunction was 
given in [9] . Obviously the above proof of shape invariance is also valid for the Hamiltonians 
of the Jl, J2 exceptional Jacobi polynomials, discussed in section [3j 

The annihilation/creation operators of the Hamiltonians of the LI, L2 exceptional La- 
guerre polynomials are obtained from those of the corresponding radial oscillator Hamilto- 
nians fl2TT6|) . fl2~T8l) . 

A^a^AKg), (5.7) 

in which are defined in (I2.8P with the replacement g — > g + £ — 1 for the LI and 
g — > g + £ + 1 for the L2. Similarly the annihilation/creation operators of the Hamiltonians 
of the Jl, J2 exceptional Jacobi polynomials are obtained from those of the corresponding 
DPT Hamiltonians (13. lip and ( I3.12p by proper replacements of the parameters. The cre- 
ation/annihilation operators for the LI exceptional Laguerre polynomials are mentioned in 
Dutta-Roy paper [17] for the £ = 1 and 2 cases. 
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The three term recurrence relations for the Laguerre and Jacobi polynomials are mapped 
to those of the exceptional orthogonal polynomials simply by Af. The explicit forms are given 
in [10]. This is definitely different from the bispectral property discussed in the previous 
section @J 
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